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Abstract 

o' 

ff ' We extend a recent work by Mussardo and Penati on integrable quantum 

field theories with a single stable particle and an infinite number of unstable 

qj \ resonance states, including the presence of a boundary. The corresponding 

scattering and reflection amplitudes are expressed in terms of Jacobian elliptic 
functions, and generalize the ones of the massive thermal Ising model and 
of the Sinh-Gordon model. In the case of the generalized Ising model we 
explicitly study the ground state energy and the one-point function of the 
thermal operator in the short-distance limit, finding an oscillating behaviour 
related to the fact that the infinite series of boundary resonances does not 
decouple from the theory even at very short-distance scales. The analysis 
of the generalized Sinh-Gordon model with boundary reveals an interesting 
constraint on the analytic structure of the reflection amplitude. The roaming 
limit procedure which leads to the Ising model, in fact, can be consistently 
performed only if we admit that the nature of the bulk spectrum uniquely 
fixes the one of resonance states on the boundary. 



1 Introduction 

The study of resonances in relativistic models has provided important insights into 
the general problem of the dynamics of quantum field theories. Resonances are 
usually due to unstable bound states which therefore do not appear as asymptotic 
states. However, associated to them, there are mass scales which may induce inter- 
esting and unexpected phenomena. Al. Zamolodchikov, for instance, has shown in 
|]]] that a single resonance state produces a remarkable pattern of roaming Renor- 
malization Group trajectories, with the property that they pass by very closely all 
minimal unitary models of Conformal Field Theory, finally ending in the massive 
phase of the Ising model. The first model with an infinite number of resonance states 
is due to A. Zamolodchikov, in relation with a QFT characterized by a dynamical Z 4 
symmetry || (see also ||). A simpler class of models with infinitely many resonance 
states have been considered by Mussardo and Penati in BJ . The analysis performed 
in Q of the Form Factors and the correlation functions of these theories has shown 
a rather rich and still unexplored behaviour of such theories. A further step in the 
analysis of infinite resonance models has been taken by Castro- Alvaredo and Fring 
in§. 

The aim of this paper is to extend the analysis of the physical effects induced 
by an infinite tower of resonance states, this time also placed at the boundary of a 
1 + 1 dimensional quantum field theory. We will consider, in particular, the boundary 
dynamics of two models which were previously defined in Q. Those may be regarded 
as the generalizations of the massive thermal Ising model and of the Sinh-Gordon 
model, respectively. In the bulk both theories have a single stable excitation of mass 
M. Our main purpose is to analyse the boundary behaviour of these theories and 
to show some novel phenomena induced by the resonance states. In particular, for 
the generalization of the Ising model we study the ground state energy and the one- 
point correlator of the thermal operator in the short-distance limit. Contrary to 
the ordinary case, where these functions scale as a constant and as an inverse power 
respectively, here they both assume an oscillating behaviour, hinting to the fact that 
the infinite series of boundary resonances does not decouple from the theory even 
at very short distance scales. The oscillating nature of these quantities prevents 
a standard conformal field theory interpretation of the theory at short distances. 
Hence, the boundary seems to have a more drastic effect on the theory with respect 
to the bulk, where the presence of an infinite number of resonances has the effect to 
make softer its ultraviolet properties leaving, however, the possibility to recover a 
conformal behaviour and to define conformal data such as the central charge or the 
anomalous dimensions 01. 
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The scattering and reflection amplitudes of the two theories investigated in this 
paper can be expressed, in the rapidity variable 9, in terms of Jacobian elliptic 
functions, which are periodic along both the complex axis. The resonance states 
correspond to poles in the unphysical sheet of the scattering amplitude. The choice 
of doubly periodic scattering amplitudes forces us to consider only theories without 
any genuine bound states, otherwise the presence of a pole on the imaginary axis 
inside the physical strip would imply a proliferation of infinite other poles located 
at the same imaginary position but at non vanishing real one, spoiling the causality 
of the theory. 

A complete list of properties of the Jacobian elliptic functions can be found in 
||. Here we simply recall that these functions, indicated by sn(it), cn(u) and dn(w), 
depend on a parameter £ (called modulus) which varies between and 1, and have 
poles and zeroes located with double periodicity on the complex plane of u. In our 
notation the rapidity variable 9 will enter the argument of these functions in the 
form u = 1^ 9, where K is the complete elliptic integral of modulus £, defined as 
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— £ 2 sin 2 a 

We will also write many expressions in terms of the complementary modulus £' = 
\J\ — £ 2 and of the corresponding complete elliptic integral K'(£) = K(£'). In the 
limit £ —* only the periodicity along the imaginary axis in 9 survives, and we 
recover the ordinary scattering theories without resonance states. The corresponding 
asymptotics of our functions are given by 

K -> - K' -> oo 
2 

sn(w) —* sin(u) cn(u) — > cos(w) dn(u) — > 1 

We anticipate here that all the physical quantities we have studied share the same 
qualitative behaviour for every £ > 0, while possible discontinous changes may occur 
for £ -> 0. 



2 Integrable Boundary Quantum Field Theories 

Integrable quantum field theories with boundary (i.e. defined on half-line and with 
the boundary placed at x — 0) have been defined and analyzed in [0, ||. Here, we 
will focus on systems with a single bulk excitation and without any bulk or boundary 
bound states. In particular, we recall that the physical strip is given by < Im 6 < i\ 
for the bulk scattering matrix S and by < lm.6 < | for the boundary reflection 
matrix K. 

In our simple case, the two basic equations of unitarity and crossing, respectively, 
assume the form: 

K(6)K(-6) = l ; (2.1) 

K(6)K(0 + m) = S{26) . (2.2) 

Following Q , where the 5-matrix was assumed to have both an imaginary and a real 
periodicity, we will look for solutions to eq. ( pTT|) and ( |2.2j) with the same property, 
i.e. 



S(6 + 2m) = S(6) 


; S(9 + T s ) = S(8) ; 


(2.3) 


<(6 + 2m) = K(6) ; 


; K(6 + T K ) = K(6) . 


(2.4) 



Eq. (|2.2|) uniquely fixes the imaginary periodicity of K to be 2iri, while allows the 
real one to be Tk = Ts or Tk = \Ts-, thanks to the fact that S is evaluated at 29; 
this observation will be useful in the following considerations. 

3 The Elliptic Ising Model 

The simplest elliptic S- matrix is given by S(6) = — 1, which satisfies eq. ( ]2.3| ) 
with any choice of Ts- Referring to Q, we interpret this as a particular analytic 
continuation of an S-matrix which possesses an infinite number of poles and zeroes, 
in the limit in which they all cancel each other. 

A solution to eq. ( |2.1|) and ( |2.2|) in correspondence to S(9) = — 1 is given by 

Km ^'Mrm- (3 ' 1) 

This elliptic function has real period T = 2n^r with the analytic structure shown 
in Figure [I], where circles and crosses represent, respectively, the positions of zeroes 
and poles. 



The poles, located at 



TX 

9 nm = — i — \- i 2rrnr + nT . 



(3.2) 



correspond to an infinite set of boundary resonance states with masses and decay 
widths given by 



M res (n) = v^M cosh (^- j , 



(nT\ 
— ) n = l,2,... 

where M is the bulk mass parameter of the theory. 



(3.3) 
(3.4) 
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Figure 1: Analytic structure of Kj(8) in the fundamental domain 



In the ordinary limit £ — > 0, Kj(9) in ( |3.1| ) reduces to the amplitude infixed (^) — 
itanh (i^ — |) found in 0. This is the reflection amplitude relative to the Ising 
model with fixed boundary conditions at x — 0. The other ordinary solution 
-Kf ree (#) = — i coth. {ij — |), relative to the Ising model with free boundary con- 
dition at x = 0, cannot be extended to the elliptic case, because it would have poles 
with real part in the physical strip. 

The elliptic generalization of the factor interpolating between Kfi xc d and K{ vee 
given in is 



sn(2Kc) - sn (if 26) sn [g(0 + lire)] en [ 



K, 



mc 



dn [f(6-iTrc 



(3.5) 



sn(2Kc) + sn (f 29) sn [g(0 - zvrc)] en [|(0 + zvrc)] dn [g(0 + zvrc) 

where the parameter c is related to the "boundary magnetic field " h by the relation 

h 2 



sn 2Kc 



2M 



(3.6) 



In order to describe the variation of h from to oo, the parameter c has to follow 
the path drawn in Figure §, which displays the analytic structure of sn(2Kc). 
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Figure 2: Path of c 



The important qualitative difference with the ordinary case is that now we have 
to impose a minimum boundary magnetic field h min = V2M, in order to start from 
the point c = in the path of Figure ^|. In fact, for c e [0, |] the factor ( |3.5|) 
has poles with real part in the physical strip, and in particular the value c = | 
(or equivalently h — 0) would correspond to the ordinary reflection amplitude Kf rcc . 
Hence, switching on the elliptic parameter I from zero to a finite value, the minimum 
admitted value of h changes discontinuously from to V2M. 

The general reflection amplitude assumes then the form 



sn r?£ (9-1%)} 

^ w= * Affirm 



IiL 

2M 



- sn (if 26) 
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2M 



+ sn(f2fl) 



^M</i<oo, (3.7) 



and in the limit h —> cxd it reduces to (|3~T 



As it was noticed, Kj(6) given by ( |3.1| ) has real period T = 27r^, but eq. (p.2|) 
with 5* (6*) = —1 admits in principle solutions with any periodicity. In particular, it is 
easy to see that the same analytic structure of ( |3~T1) , but with half period T" 
is realized by 



7T 



K ' 



K'M 



^ifA^-m 



m [K(6l-<|)]dn[£(e-if)] • 
Since it will be used in the following, we also note that the function 

K'l(6) = -zK' I (e)K' I (e--' 



(3.* 



(3.9) 



has again the same structure, but with real period T" = |^ . 

It is worth to stress, however, that the various choices for the period of the 
reflection amplitude have a different physical meaning, corresponding to theories 
with distinct spectra of boundary excitations (see eq. ( |3.3Q and ( |3.4j) ). 

3.1 Ground State Energy 

In || it was developed the technique to compute the ground state energy in in- 
tegrable systems with boundary. In particular, in the case with S(9) = —1 the 
resulting expression for the effective central charge with equal boundary conditions 
at both sides of the strip takes the simple form: 

12 f°° / \ 

c cS = — limr/ d0cosh01n(l + |A:(0)| 2 e- 2rcoshO ) , (3.10) 

where K{6) = K (Jf — 9j is the amplitude relative to the boundary condition under 
consideration. We have studied the behaviour of this quantity choosing the elliptic 
reflection amplitude^] Ki{9) in Q3.1| ). 

Performing the change of variable x = r cosh 9 in order to calculate the integral 
in ( p.lO| ), it is easy to see that the limit r — ► implies an evaluation of \Kj(9)\ 2 at 
9 — > oo. In the ordinary case this simply corresponds to \K &xcd (9)\ 2 —> 1, but in the 
elliptic case \Ki(9)\ 2 is a periodic oscillating function also for real values of 9, and 
this has the consequence that the limit procedure does not converge to a fixed value 
of c c fj. In particular, if we evaluate the limit r — > on the sequence {r„} n=1 2 , 

with r n = TTTrX i the corresponding sequence of integrals converges very rapidly to 

some fixed value which, however, depends on the choice of a. This indicates that the 
scaling function c e s(Mr), which in the ordinary case tends to a constant as Mr — > 0, 
keeps a dependence on the variable Mr, assuming the same values for constant 
values of M res (n)r n . This situation has a clear physical interpretation in terms of 
the presence of boundary resonance states with infinitely increasing energies, which 
do not decouple from the theory even at very high energies. 

Figures § and f| show the behaviour of c c s as a function of the parameter a for 
two values of I, i.e. £ = 0.32 and £ = 0.95, respectively. In both cases, the first 
plot corresponds to a variation of a from to 10, while the second one describes the 
variation of w between 2 and 10, where we have defined a = 10 -u \ The logarithmic 
scale clearly displays the oscillating behaviour, which leads to the absence of a given 
limit. 



1 Alternatively, we could have chosen Kj(9) or Kj(9) given, respectively, by (3.8) and (3.9), but 
for convenience, we have preferred to use the simplest expression. 





Figure 3: I 
w G [2, 10]. 



0.32. (a) c c fj vs a, with a G [0, 10], (b) c e fj vs w, with a = 10 w and 





Figure 4: £ 
w G [2, 10]. 



0.95. (a) c c ff vs a, with a G [0, 10], (6) c c ff vs w, with a = 10 w and 



3.2 One-point Function of the Energy Operator 

It is equally interesting to analyse the behaviour of the correlators in the presence of 
an infinite number of boundary resonance states. Here we will consider the simplest 
correlation function. 

The expression for the one-point function of the energy operator in the thermal 



Ising model with boundary has been studied in |L0 |, and it is given by 



e (t) = ^2<0\e(x,t)\n><n\B >, (3.11) 

n=0 

where | n > is the asymptotic bulk state with n particles, and the boundary state 
| B > is given by a superposition of two-particle states in which every couple consists 
of excitations with opposite rapidities: 

ro ° dB 



B >= exp 



uo 2tt 



K{0)A\-0)A\0) 



> . (3.12) 



Using the corresponding elliptic form factor calculated in ||, we get the following 
expression for the matrix element of the energy operator on the two-particle state, 
which is the only non vanishing contribution: 



<O|e(s,t)|0i,0 2 > = 2ttM^ ^%^ 6 *\\ : 

dn 1.^(01-02 J 



exp [-M£(cosh0i + cosh0 2 ) + iMx(sinh0i + sinh0 2 )] . 

The correlator assumes then the form 

poo sn (^20) A 
e (t) = I'M / d6 — W—v K (0) e - 2Micosh9 . (3.14) 

As before, we choose to perform the calculation with the reflection amplitude Kj[ff) 
given by ( 3-l| ); hence, eq. (|3.14j ) specializes to 

e (t) = -M(£') 3/2 / d6E(e)e~ 2Mtcoshe , (3.15) 

Jo 



where 



, s sn (£20) sn (£0) 

BW = d^t4'^fcy- (3 - 16) 



We are interested to study ( |3.15| ) in the short- distance limit. Contrary to the ordi- 



nary case (£ — 0), where this one-point function scales as a power law JlQfl , now we 
expect at most a logarithmic divergence, due to the fact that E(6), being a periodic 
function without poles on the real axis, is limited on the range of integration. 
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However, as it is shown in Figures ^| and ^ (with i = 0.71), the correlator neither 
diverges nor goes to a finite limit as Mi — > 0, but it has an oscillating behaviour, 
more visible on the logarithmic scale. 




Figure 5: e (r) for r G [0, 1] (r = 2Mt), 



0.71 




Figure 6: e (w) with r = 10 _u \ w G [2, 10], 



0.71 



This can be easily understood analyzing the interplay between the two different 
factors present in the integrand. In fact, E(6) oscillates assuming positive and 
negative values with the same amplitude (Figure [7]), while for small values of r the 
exponential exp(— r cosh^) is approximately constant on an interval in 9 of order 
In -, and then falls rapidly to zero. Hence, if we decompose the length of this interval 
as In - = nT + x, the integration on nT gives roughly zero, while the remaining part 
corresponds to an integration of E{6) on a certain fraction of its period, so that it 
oscillates depending on the value of x. 




Figure 7: E(8) with £ = 0.71 

Let us now turn to the reflection amplitude (|3.7| ) with the boundary parameter 
h. In this case the one-point function is given by ( p.!5|) with 



E{9) 



l-^-«n[^(»i-6>)] 
sn [2? (if - 0)] 



IiL 

2M 



E{9). 



(3.17) 



For every h m [ n < h < oo, at short distances these correlators do not tend to a 
given limit, as before. Figure |8| shows the one-point correlator as a function of r 
for different choices of h, in the case £ = 0.55. We note that the value h = h min 
corresponds just to a global change of sign with respect to case of Ki{9) (equivalent 
to h — > oo). 

Summarizing, the common behaviour shared by all these correlators is not sur- 
prising in the light of the previous result about the dependence of the central charge 
on the scaling variable Mr, even at short distances. Also in this case, the absence 
of a proper ultraviolet limit can be related to the fact that the infinite boundary 
resonance states never decouple from the theory. By using the Form Factors of the 
magnetic operator computed in M , it is also easy to see that the same phenomenon 
occurs for its one-point function of this field. 



10 



-0.25 

-0.5 

-0.75 

-1 

-1.25 

-1.5 
-1.75 



0.2 0.4 0.6 0. 



(b) 



0.0025 




(d) 




0.05 



0.2 0.4 0.6 0.8 1 



Figure 8: £ = 0.55. One-point function of the energy operator with (a) h = h r 
(b) h = 1.60 /i m in, (c) h = 1.61 hmia and (d) h = 1.73 h min . 
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4 The Sinh-Gordon Model 

4.1 The Ordinary Case 

The Sinh-Gordon model is defined by the field equation 



dU - dlcf) + y^f- sinMVW) = , (4.1] 



and has a bulk S'-matrix of the form ( [[LI]] ) 



where 



s a {0)si„ a (6) 



= smh[l(6 + nrx)] ^ 
sinh [t}(0 — Z7rx)l 



and 

(3 2 



(4.4) 



4vr Vl + /9 2 /4vr / 

In the presence of a boundary, this model is integrable if the boundary conditions 
have the form (0,0) 



dA\, = ^f [e e-^^> - e^**™) , (4.5) 

where Eq and £i are two boundary parameters, which can be alternatively parame- 
terized by 

Si = cos ircii < ai < ao < 1 . (4-6) 

The reflection amplitude can be obtained by analytic continuation from the one 
of the sine-Gordon breather calculated in ||13| , and has been analyzed in |TJ]]. It is 
given by 

si(6)si+a(e)s2- a (9) 
K(6) = 2 2 2 (4 7) 

where 

E= (a + ai)(l - a) ; F = (a - ai)(l - a) . (4.8) 

In the case of Neumann boundary conditions, where ao = a,\ = ~, eq. ( |4.7| ) becomes 

Sl±a(0) 
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4.2 The Elliptic Case 

Let us consider the elliptic S'-matrix^ 

L in \ 



S{6,a) 



sn 



ma en 



Lin ' 



27ra 



sn[£(# + ZTra)]cn[£(fl-z7ra 



(4.10) 



which has real period T = 27r^ and in the ordinary limit (£ —> 0) reduces to the 
Sinh-Gordon S'-matrix ( [1-.2| ). 

The analytic structure of S(#, a) in shown in Figure §. 
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Figure 9: Analytic structure of S(9) in the fundamental domain 



In Q] it was analyzed the generalization of the roaming limit, which leads to the 
S'-matrix of the Ising model S(8) = — 1. In particular, it consists in taking 



T 
a = m + i(2n + 1) h S 

2n 



m G Z, n G Z 



and sending 5^0. The arrows in Figure ^| refer to the case m = n 
the final positions of poles and zeroes, which cancel each others. 



(4.H) 



and indicate 



2 In H it was constructed and studied an elliptic S'-matrix with a very similar structure, but 
with half real period T" = ff-fr- Our present choice is simply for convenient reasons. In fact, by 
fixing the imaginary period of S and K to be 27ri, the only two real periods which can be easily 
obtained multiplying Jacobian elliptic functions are exactly T and T'\ a period T", for example, 
can be obtained only in the more subtle way related to ( |3.9| ). As we will see below, the K (0) is 
forced to have half the real period of the S'-matrix. Hence, starting with the S matrix of H , we 
would have had from the beginning a cumbersome expression for the reflection amplitude, which 
we prefer to avoid. 
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We present now an elliptic reflection amplitude with real period X" = 7r|^ and 
without poles having a real part in the physical strip, which in the ordinary limit 
tends to the Neumann reflection amplitude of the Sinh-Gordon model (|4.9|): 



K (9) = K'}{9)f{9)N{9) 
where K"{6) is given by ( |3.9|) , the factor 

«a[g(*-*fa)Mg (* + *!(! 



N{9) 



0)] 



•1(1 



is the direct generalization of the ordinary one, and 

^[l i (9-qa)]^[l l (9 + q(l + a))} 



))] 



W) 



^[i i {9 + q a )]^[j- i {9-q(i + a))} 



(4.12) 



(4.13) 



(4.14) 



is a doubly periodic function necessary to satisfy the crossing equation (|2.2j ), which 
tends to 1 as £ — > (for the properties of Theta functions, see 0). 
The analytic structure of K{9) is shown in Figure [l(^ 
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Figure 10: Analytic structure of K{9) in the fundamental domain 



It is interesting to note that the initial Ising factor present in (|4.12j) is the one 
with period T" = f ^-, while ( f4.12| ) has altogether period T' = vr^. This choice, 
which may not seem very natural, is motivated by the requirement that the roaming 
limit applied on (|4.12|) produces the reflection amplitude of the Ising model. 
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In fact, considering a = m + i(2n + l)j--\- 5, in the limit § —> poles and zeroes 
move according to Figure |ll|, leading to the Ising reflection amplitude (|3.8|) with 
period T'. 
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Figure 11: Roaming limit ( 



O 



-T' 



x 



- > *£) 



ITT 



O 



— 111 



:< 



[£ 



o 



T' 



x 



It is now clear the reason why we have to use a reflection amplitude with real 
period half of the S'-matrix one. In fact, contrary to the ordinary roaming limit in 
which a = 1 + iC and C — > oo, now the imaginary part of a is sent to a finite value 
2~. Hence, since in the reflection amplitude the parameter a is multiplied by a factor 
| absent in the expression of the S'-matrix, poles and zeroes are shifted on the real 
axis by ^ instead of -|, and, in order to reach the same position, at the beginning 
they have to be located at half the initial distance. 

Figure [H] refers to the case m = n = 0, but the same final situation will take 
place for any n G Z and for m = 0,3 (mod 4). On the other hand, in the case 
m = 1,2 (mod 4) we need to consider an initial K(9) with an Ising factor Kj(9) 
instead of K"(9), recovering at the end K'({9). 



In principle, it is possible to extend to the elliptic case also the general reflection 
amplitude (|4.7|) , defining 



K{9) = K'l(9)f(9)G(9) 



(4.15) 
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where 



G(0) 



m [%(e + q(i + a))]m[%(e + q(2-a))] 



sn 



[g(e-<|(l + a))] m [g(e-ii(2-a))] 

^[SP-»'f(i--s))]™[SP-»f( 1 + ^))] 



^[f (^-iKl-F^sng^-iKl + F))] 



^[gp + «'i(i-^))]^[gp 



j|0_ 
<f(l 



*"))] 



(4.16) 



and 
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i(G> + if(l + a))]04[£(0 + ti(2-a))] 
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■»f(l + 0))]^4 
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[^(e + qii + E))] 


#4 


[£('- 


■tf(l- 


- F))] tf 4 


[£('- 


-tf(l + F))] 



l(0 + if(l-F))]tf 4 [i(0 + if(lH-F))] 



(4.17) 



In order to avoid poles with real part in the physical strip, both the parameters 
E and F have to lie in the interval [0, 1]. Furthermore, the roaming limit leads to 
an Ising reflection factor only for certain specific choices of E and F in terms of a. 



5 Conclusions 

In this paper we have studied the physical effects induced by an infinite number of 
boundary resonance states in certain integrable quantum field theories. The main 
result of our analysis regards the ultraviolet regime shown by such kind of theories. 
We have studied, in particular, in the elliptic generalization of the Ising model, the 
short-distance behaviour of the effective central charge and of the one-point function 
of the energy operator, which in the ordinary case scale, respectively, as a constant 
and as -. In the elliptic case, instead, both these quantities display oscillations as 
r — > 0, and do not have a definite limit. This leads to the physical interpretation 
that the infinite resonance states with increasing mass present on the boundary 
never decouple from the theory, even at very short distances, making meaningless 
in this case the concept of scaling limit. 

The other interesting feature emerged from this work is the crucial role played 
by the choice of the real periodicity. In fact, various theories with different real 
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periods for the scattering and the reflection amplitudes lead to the same theory 
in the limit 0. — > 0, because the ordinary case corresponds to T — > oo. However, 
since for £ > this period is finite, it is necessary to treat it carefully. Indeed, the 
roaming limit procedure, which connects the Sinh-Gordon model to the Ising model, 
can be consistently generalized to the elliptic case only imposing that the reflection 
amplitude has half real period with respect to the S'-matrix. This has a physical 
implication on the nature of the spectra of the theory: once the bulk spectrum of 
resonance states is given, the boundary one is then automatically fixed. 
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